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ABSTRACT 

C3 . We present multi-colour photometry of the M8.5V ultracool dwarf "pulsar" TVLM 

513-46546 (hereafter TVLM 513) obtained with the triple-beam photometer ULTRA- 
I , CAM. Data were obtained simultaneously in the Sloan-g' and Sloan-i' bands. The 

■ previously reported sinusoidal variability, with a period of 2-hrs, is recovered here. 
£f} ' However, the Sloan-g' and Sloan-i' lightcurves are anti-correlated, a fact which is in- 
0^ . compatible with the currently proposed starspot explanation for the optical variability. 
t-H ■ The anti-correlated nature and relative amplitudes of the optical lightcurves are con- 

| sistent with the effects of persistent dust clouds rotating on the surface of the star. In 

, the absence of other plausible explanations for the optical variability of TVLM 513, it 

' seems likely that dust cloud coverage combined with the rapid rotation of TVLM 513 

OO ! is responsible for the optical variability in this object. However, crude modelling of a 

■ photosphere with partial dust cloud coverage shows that the anti-correlation can only 

be reproduced using cooler models than the literature temperature of TVLM 513. We 

i-h , suggest this discrepancy can be removed if more dust is present within the photosphere 

>^ ■ of TVLM 513 than theoretical model atmospheres predict, though a definitive state- 

ment on this matter will require the development of self-consistent models of partially 
dusty atmospheres. 
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1 INTRODUCTION 

Brown dwarfs and very-low mass stars (collectively known as 
ultracool dwarfs, or UCDs) are strongly affected by the pres- 
ence of dust in their photospheres. Dust absorbs elements 
from the gas phase, changing the opacity and metallicity of 
the atmosphere. Dust begins to form at temperatures corre- 
sponding to the transition between M and L spectral types, 
and becomes more prominent as the atmosphere cools; the 
presence of dust thus defines the L-dwarfs. Finally, when the 
dust clouds "rain-out" at lower temperatures still, they are 
responsible for the L-T transition. Thus understanding the 



1 Based on observations made at the European Southern Obser- 
vatory, Paranal, Chile (ESQ program 079.C-0686) 



formation, chemistry and atmospheric dynamics of dust is 
the central challenge facing the ories of ultracool dwarf at- 
mospheres |Burrows et al. 2006h . 



The presence of dust is also thought to affect the mag- 
netic properties of ultracool dwarfs. A combination of an 
increasingly neutral atmosphere, and frequent collisions be- 
tween charged particles and dust grains in the dense at- 
mosphere results in the atmosphere s of ultracool dwarf s 
having a high electrical resistivity (|Mohantv et al. 20021 ). 
Thus, whilst strong quiescent and flaring radio emis- 
sion reveals strong magnetic fi elds amongst the ultracool 
dwarfs (e.g Berger et al. 200ll : iBurgasser fc Putman 20051 ; 
lHallinan et al. 20071 ). the predominantly neutral atmo- 
spheres may explain the relatively low levels of other ac- 
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tivity indicators , such as Ha em ission and X-rays (e.g 
iGizis et al. 200d : IWest et al. 20041 ). 

TVLM 513 (M8.5V) is an ideal object in which to 
study the interplay of magnetic fields wi th a cool, dense at - 
mosphere. It is close by (d = 10.6 pc; iDahn et al. 2 002). 
has known radio and Ha activity, and is rapidly rotat- 
ing (|Mohantv fc Basri 200$ ). The radio emission of TVLM 
513 is fascinating; observations from 1.4 to 8.5 GHz show 
the radio emission can switch st ates between highly po- 
larised pulses with a 2-hr period l|Hallinan et al. 20071 ) to 
a fainter, quiescen t state interrupted by stochastic flares 
l|Berger et al. 20081 ). The periodic nature of the radio emis- 
sion has led to this object being dubbed an ultracool 
dwarf "pulsar". Strong Ha emission, modulated on the 
same 2-hr period , and a detection of X-ray emission 
l|Berger et al. 20081 ). show that TVLM 513 is a magnetically 
active star, which is capable of supporting a chromosphere 
and a coron a. This is despite its low effective temperature 
of -2300 K (|Dahn et al. 20021 ). roughly the temperature at 
which dust is believed to start forming within the photo- 
sphere fsee lBurrows et al. 20061 . and references within). 

Optical I-band photometry of TVLM 513 revealed the 
same 2-hr periodicity as seen in the radio and Ha emission ; 
this was assumed to be due to starspots (|Lane et al. 2007T ). 
However, optical variability in ultracool dwarfs can also 
be caused by du st clouds within the photosphere (see 
iBailer- Jones 20021 , for a review). Simultaneous optical pho- 
tometry in two or more bands can distinguish between vari- 
ability cause d by starspots (|Rockenfeller et al. 20061 ) and 
dust clouds (jLittlefair et al. 20061 ). We therefore obtained 
simultaneous, multi- colour photometry using the fast pho- 
tometer ultracam (|Dhillon et al. 20071 ) to determine the 
cause of the optical variability in TVLM 513. The obser- 
vations are described in Section [5] the results presented in 
Section [3] and discussed in Section [4] whilst in Section [5] we 
draw our conclusions. 



2 OBSERVATIONS & DATA REDUCTION 

For a duration of 2.7 hours on the night of June 9 th 2007, the 
M8.5V star TVLM 513 was observed simultaneously in the 
Sloan-p' and i' bands using ultracam on the 8.2-m melipal 
unit of the Very Large Telescope (VLT) at Paranal, Chile. 
Thanks to the frame-transfer CCDs employed in ultracam 
we were able to obtain ~10 second exposures with 24 msecs 
dead-time between frames. Sloan-u' images were obtained si- 
multaneously with the other colours, but the star is too faint 
to produce useful lightcurves in this band. The Sloan-g' and 
i' magnitudes of TVLM 513 are 20.9 and 16.1, respectively. 
Observing conditions were excellent; the night was photo- 
metric and the seeing stable around 0.6 arcseconds. Data 
reduction was carried out using the ultracam pipeline data 
reduction software. The data were first bias subtracted and 
flat fielded. The flat field used was a median of ~100 twi- 
light frames taken on a blank piece of sky with the telescope 
dithering. Because the flat fields suffer from scattered light, 
large-scale trends were removed from the flat-field before 
use, by dividing the flat-field by a median-filtered version of 
itself. Since the position of the stars on the CCD remains 
constant throughout our observations, this will not affect 
our photometry. 



Extraction of target and reference star lightcurves was 
performed using a variable-sized aperture. The aperture size 
for both objects was set to 1.7 times the full-width at half 
maximum of the target stars spatial profile, as measured by 
a Gaussian fit. Differential photometry was obtained by di- 
viding by the flux of a nearby reference star. This implicitly 
assumes that the same proportion of the total flux falls out- 
side the aperture for both target and reference stars. This 
is only true if the target and reference stars have the same 
point spread function (PSF). If the instrumental PSF of the 
target and reference stars differ, this introduces errors into 
the photometry. Worse still, since the actual PSF is the con- 
volution of the instrumental PSF and the seeing profile, it is 
true that the difference in PSF between the target and ref- 
erence (and hence error in the photometry) is a systematic 
function of seeing. To minimise this source of error one can 
either use many reference stars, and fit the PSF as a function 
of chip position and time, or you can use a single reference 
star as close to the target as possible. Our images do not 
contain enough reference stars to implement the former ap- 
proach, and so here we choose to implement the latter. The 
reference star chosen has g' and i' magnitudes of 17.2 and 
16.6 and is 37 arcseconds from the target star. 

Photometry for multiple reference stars was extracted 
and we have used this both as a "sanity check" and as a mea- 
sure of our photometric accuracy (although, for the reasons 
discussed above, this is likely to be an upper-limit). The pho- 
tometric accuracy we determine is 3 mmag for the i'-band 
and 7 mmag for the (/-band. This is smaller than the sinu- 
soidal variability reported in section [3] giving us confidence 
that this variability is real. Indeed, the sinusoidal variability 
is recovered, albeit with more noise, if other reference stars 
are used to produce the differential photometry. 



3 RESULTS 

The ultracam photometry is shown in Figure [T] Both the 
g' and i' lightcurves exhibit sinusoidal variability with a 
period of around two hours. To assess the significance of 
this variability, 1000 copies of the data were made with 
the time-ordering of the data points randomised. Lomb- 
Scargle periodograms of these copies were computed and 
compared to the periodograms of the original data. None of 
the randomised datasets showed power in the periodograms 
as large as the peak power of the actual data, and so the false 
alarm probability for the periodic variability is less than 0.1 
per cent in both bands. Errors in the periods were com- 
puted using a bootstrapping technique. We computed 1000 
copies of the data using bootstrapping with replacement. 
Periods for these copies were determined from the peak of 
a Lomb-Scargle periodogram, and the scatter in these val- 
ues used as an estimate of our period error. We find peri- 
ods of 2.0 ± 0.5 hours for the (/-band and 2.4 ± 0.4 hours 
for the i'-band. This confirms the perio dic variability see n 
in the I-band lightcurves of TVLM 513 l|Lane et al. 20071 ). 
which is also seen in the equivalent width of Ha emission 
|Berger et al. 20081 ). 

The RMS variability in the (/-band is about 1.5 per 
cent. RMS variability in the i'-band is 0.15 per cent. The i'- 
band variabilit y is thus much sm aller than that reported in 
the I-band bv lLane et al. f2007l ). which was nearer 0.8 per- 
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Figure 1. Sloan-c/ and i' lightcurves of TVLM 513. The change in brightness of the target star (mmags), is plotted against time in 
days. The data have been binned into bins of 20 second duration. The Sloan- g' lightcurve is plotted in green in the top panel, whilst the 
bottom panel shows the Sloan-i' lightcurve, plotted in red. Note the vertical scales for the (/-band lightcurve is much larger than that 
for the i'-band lightcurve. Also shown are sinusoidal fits to the data 



cent. The most striking fact about the optical variability is 
that the g lightcurve and i lightcu rves are in anti-phase. A 
discrete cross-correlation function l|Edelson fc Krolik 19881 ) 
yields a phase shift between the two lightcurves of 0.50±0.05. 
It is not easy to reconcile this observation with the claim 
that the bro adband optical va riability of TVLM 513 is due 
to starspots l|Lane et al. 20071 ); since starspots cause a dim- 
ming of the star at all wavelengths, starspot-induced vari- 
ability would produce g' and i' lightcurves which are in 
phase. We therefore turn our attention to other possible 
causes of the optical variability. 

3.1 Starspots and Faculae 

Whilst starspots alone cannot reproduce the optical variabil- 
ity, it is possible in principle to explain it using a combina- 
tion of starspots and faculae. Because the faculae are hotter 
than the surrounding photosphere, we can conceive of a sit- 
uation in which the effect of cool starspots dominates in the 
i'-band, whilst the faculae have a dominant effect in the g'- 
band. In these circumstances, an excess of spots would cause 
a dimming in the i'-band. Due to the corresponding excess 
of faculae, this would be accompanied by a brightening in 
the (/-band. 

In the absence of constraints on spot and faculae prop- 
erties in ultracool dwarfs, we cannot rule out this scenario 
entirely. However, we note the following points, which make 
this explanation unlikely. First, to ensure that the effect of 
spots dominates in the i'-band, whilst the effects of faculae 
dominate in the (/-band requires considerable fine-tuning of 
the relative areas of spot and faculae coverage. We find, that 



for plausible spot and faculae temperatures of T s — 1800 K 
and Tf = 2400 K, respectively, the ratio of faculae area to 
spot area must lie between 1.4 and 2.2. This fine tuning 
of relative areas makes this scenario unlikely a-priori. Fur- 
thermore, the g' and i' lightcurves are observed to vary in 
anti-phase to a good level of accuracy (see Figure [1} • This 
requires additional fine-tuning of the spatial arrangement of 
the spots and faculae, as the contrast of faculae is greatest 
at the limb, wh ereas spot cont rast is roughly independent 
of limb position |Gondoin 20081 ). For these reasons, we con- 
sider an explanation of the optical variability in terms of 
spots and faculae unlikely. 

3.2 Chromospheric emission 

Anti-correlated variability in the Sloan-</ and i! bands is 
possible with a combination of starspots and chromospheric 
emission. Whilst the CaHK lines lie outside the g' and i' 
bands, several of the higher-order Balmer lines fall in the 
g' band. Thus a more active region of the star might show 
lower i'-band flux (due to the dimming caused by starspots) , 
but a higher (/-band flux (due to the contribution to the g'- 
band from Balmer emission lines). However, the spectrum 
of TVLM 513 appears to rul e out this conjecture (Figure^). 
From the Gemini spectra of iBerger et al. (20081 ) we can see 
that the higher-order Balmer lines contribute very little to 
the (/-band flux, even in the high state. In fact, Balmer 
emission only contributes an additional 0.3 per cent to the 
(/-flux, so it is unable to account for the 3 per cent peak- 
to-peak variability seen i n the g'-lightcurve. Since the Ha 
equivalent widths in the IBerger et al. f200ct ) spectra span 
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Figu re 2. Ge mini spectra of TVLM 513 from 
iBerger et al. (200t$) . The spectra in high and low-Balmer 
emission states are shown. The high Balmer emission state is 
plotted in black, with the low-Balmer emission state shown in 
light grey. Both spectra have been normalised by the peak flux 
of the high Balmer emission state spectrum. The ULTRACAM 
response curve in the g'-band is overplotted. A weak H/3 emission 
line is visible in the high-state spectrum. 



the range of observed Ha equivalent widths in the literature 
for TVLM 513, it is likely that the Balmer emission in these 
spectra adequately reflects the Balmer emission at the time 
of our observations, and so we judge it unlikely that chro- 
mospheric emission is responsible for the optical variability. 

3.3 Magnetospheric Emission 

The radio emission of TVLM 513 is highly variable in nature. 
The radio emission seems to change b etween the bright pul - 
sations with 2-hr period reported bv lHallinan et al. (20071 ) 
to a fainter, quiesce nt, state interrupted by stochastic flares 
|Berger et al. 20081 ). There is some debate as to whether 
the radio emission results fro m an electron cyclotro n maser 
(ECM) emission mechanism (Ha llinan et al. 20061 ). or the 
same gyrosynchrotron emi ssion which powers the radio emis- 
sion in warmer M-dwarfs (jBerger et al. 20081 ). As a narrow- 
band emission mechanism, ECM emission will not con- 
tribute to the optical flux, and so here make the conservative 
assumption that the radio emission is caused by gyrosyn- 
chroton radiation. 

It is beyond the scope of this letter to determine 
whether gyrosynchrotron radiation can provide sufficient 
flux at optical wavelengths to cause the variability seen here. 
However, since gyrosynchrotron radiation produces a broad 
continuum, it is not likely that it is directly responsible for 
the anti-correlation between the g' and i' band variability. It 
is possible to construct a viable physical model whereby the 
i'-band variability is caused by starspots, and the (/-band 
variability is caused by variations in the gyrosynchrotron ra- 
diation. Such a model can plausibly produce anti-correlated 
variability between the two bands, if areas of the star with 
an excess of starspots are associated with brighter gyrosyn- 
chrotron radiation. The sinusoidal nature of the (/-band 
lightcurve argues against such a model, however, as the emis- 
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Figure 3. Top: Calculations of the change in spectrum of an ul- 
tracool dwarf, due to the appearance of a dust cloud that covers 
10 per cent of th e stellar surface. T he calculations use the syn- 
thetic spectra of lAllard et al. f200lh - see text for details. The 
thick solid line shows a model with T e g = 1900 K, whilst the 
thick dashed line shows a model with T c ff = 2300 K. The thin 
solid line shows a hybrid model in which dust forms within a 
1900 K starspot on a 2300 K photosphere. Bottom: Also shown 
are the bandpasses of the Sloan-g' filter and the Sloan-i' filter. 



sion in the radio is not observed to vary sinusoidally. We 
therefore believe it is unlikely that magnetospheric emission, 
or a combination of magnetospheric emission and starspots, 
is responsible for the optical variability. 

3.4 Dust clouds 

The presence of dusts clouds can, in principle, explain the 
optical variability seen in TVLM 513. Consider the forma- 
tion of a dust cloud in a largely dust-clear photosphere. The 
cloud increases the continuum opacity, whilst at the same 
time the condensation in dust clouds of gas-phase alkali- 
metals reduces the opacity in spectral regions dominated by 
alkali absorption. Since the Sloan-p' band is dominated by 
continuum opacity, whilst the Sloan-i' band is dominated 
by molecular absorption, the formation of a dust cloud will 
reduce (/-band flux, and simultane ously increase the i'-band 
flux. Following Bail er- Jones (2002) we create a toy model of 
this p rocess by combinin g the COND and DUSTY0 mod- 
els of lAllard et al. f200ll ) in the ratio 1:9, and comparing 
the resulting spectrum with a 100 per cent COND model 
of the same effective temperature. We would like to stress 
that this model is only meant as an illustration, and is far 
too crude to attempt a detailed analysis such as lightcurve 
modelling, or derivation of cloud coverage. Such an analysis 
would require more sophisticated modelling which properly 

1 Both COND and DUSTY models include dust for the purposes 
of calculating the chemistry, but COND models assume the dust 
settles below the photosphere, whilst DUSTY models assume the 
dust remains in the visible part of the atmosphere 
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and self-consistently treated the dust content of cloudy and 
"cloud- free" regions of the photosphere. Examples of the 
spectral changes resulting from our toy model are shown in 
Figure [3] for effective temperatures of 1900 K and 2300 K, 
as compar ed to the effective temperature for TVLM 513 of 
~2300 K ijDahn et al. 20021 ). 

The models in Figure [3] show that both the g' and i' 
bands show spectral changes due to resonance lines and 
changes in the continuum brightness. Whether the star gets 
brighter or fainter in a given band depends on whether line 
or continuum changes dominate. We determine this by fold- 
ing the spectral changes in Figure [3] through the g and i' 
bandpasses. For the 2300 K model we find the continuum 
change dominates in both bands; the warmer model thus 
predicts correlated variability in the g' and i' bands. By 
contrast, in the 1900 K model, the resonance lines dominate 
the j'-band whilst the continuum dominates the (/-band. 
The cooler model thus predicts the observed anti- correlated 
g and i' variability. Therefore, only by using significantly 
cooler synthetic spectra than appropriate for TVLM 513 can 
we obtain the observed anti-correlation between g' and i' . 
Figure[3]shows why this should be the case. Whilst the mod- 
els at 2300 K and 1900 K produce spectral changes of similar 
shape, it is clear that, in the 2300 K model, the change in the 
alkali metal lines is much smaller. This is presumably a nat- 
ural consequence of the fact that the amount of dust formed 
in the warmer model is itself much smaller than in the cooler 
one. Therefore, if the models underestimate the amount of 
dust that forms in a 2300 K photosphere, it would be possi- 
ble to reproduce the observed anti-correlation variability in 
g' and i' using models of the appropriate effective tempera- 
ture. One might suspect this problem can be overcome using 
a "hybrid model" in which dust clouds form in the cooler 
regions of the photosphere created by starspots. Figure [3] 
shows that this is not the case. In the figure we show such a 
model, with a dusty, 1900 K starspot covering 10 per cent of 
the surface of a dust-free 2300 K photosphere. The starspot 
introduces such strong variability that it overwhelms the ef- 
fects of the dust clouds, and this model produces correlated 
variability in the g' and i' bands. 



variability is not dominated by starspots, as was previously 
assume d. Indeed, optical multi- band photometric monitor- 
ing by iRockenfeller et al. (20061 ) shows that starspots can 
cause periodic variability in objects as late as M9V. In this 
context, it i s worth mentionin g that the I-band variability 
reported bv lLane et al. (20 071 is four times larger than the 
variability reported here in i . Such large amplitudes are 
difficult to explain using a dust cloud model, and may in- 
deed be due to starspots. Thus, there is some evidence that 
the nature of the optical variability in TVLM 513 might 
not be stable, being at times dominated by dust clouds, 
at other times by starspots. Further multi-wavelength ob- 
servations will be needed to see if the nature of the opti- 
cal variability is linked to changes in t he radio behaviour 
l|Hallinan et al. 20071 ; iBerger et al. 20081 ) of this remarkable 
object. 



5 CONCLUSIONS 

We present simultaneous monitoring of the M8.5V star 
TVLM 513 in the Sloan-p' and i! bands. Both bands show 
sinusoidal variability on the ~2-hr rotation period. The 
g' and i band lightcurves are in anti-phase, a fact which 
is incompatible with the optical variability being due to 
starspots. Without a plausible alternative explanation for 
this behaviour, we conclude that the optical variability in 
TVLM 513 is caused by the presence of a dust cloud in a pre- 
dominantly dust-free photosphere. The optical variability is 
consistent with a dust cloud covering a significant fraction 
of the photosphere, but only if cooler stellar models than 
appropriate are used. This could be a consequence of the 
crudity of our modelling, or it may imply there is more dust 
in the at mosphere of TVLM 5 13 than the stellar atmosphere 
models of All ard et al. (200ll ) predict. The i'-band variabil- 
ity reported here is a factor o f four smaller than previously 
reported bv lLane et al. (20071 ). Since such a high amplitude 
of variability is difficult to reconcile with dust cloud models, 
we suggest that the optical variability of TVLM 513-46546 
can change in nature. It is unclear if this is related to the 
observed changes in radio emission from this object. 



4 DISCUSSION 

Based upon a lack of plausible alternative explanations for 
the optical variability in TVLM 513, we suggest that the 
most likely cause of the variability we observe is the rota- 
tion of a dust cloud covering a significant fraction of the stel- 
lar photosphere. This explanation is not without its prob- 
lems. Crude modelling of a photosphere with partial dust 
cloud coverage shows that the observed optical variability 
can only be reproduced using models which are too cool 
by so me 400 K. This ma y be an indication that the mod- 
els o f lAllard et al. (2001) underestimate the amount of dust 
that has formed in the photosphere of TVLM 513. Alterna- 
tively, this discrepancy could be a consequence of our sim- 
plistic modelling process; a definitive statement on this mat- 
ter awaits the arrival of self-consistent models of the emission 
from partially dusty atmospheres. 

Based upon the strong magnetic field inferred from ra- 
dio observations of TVLM 513, and the strong Ha emission 
seen in this object, it is perhaps surprising that the optical 
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